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NASA TECHNICAL PAPER 


SPACE SHUTTLE MAIN ENGINE HIGH PRESSURE FUEL PUMP 
AFT PLATFORM SEAL CAVITY FLOW ANALYSIS 

I. INTRODUCTION 


In working to improve the performance of the Space Shuttle Main Engine (SSME), the engineer is confron- 
ted with the difficult task of analyzing a complex engine system running under extreme operating conditions. The 
temperatures in the Shuttle’s engines range from 37°R up to 5000°R, pressure vary from 20 to 8000 psi, and the 
engines’ pumps rotate at speeds up to 37,000 rpm. Direct measurement of the engine environment is often impracti- 
cal. Indeed, the particular area under consideration may be virtually inaccessible to instrumentation. Fortunately, 
the capability of modeling heat and mass transfer using computers has advanced to the point where computational 
fluid dynamics (CFD) can provide an alternate method of analyzing the engine. When used to model the various 
components and processes in the engine, numerical analysis can provide the engineer with valuable insight by 
allowing him or her to examine a wide range of operating conditions. The effect of a change in geometry, of a 
change in flowrate, or of a change in any parameter can be examined. Even a simple numerical model can demon- 
strate the sensitivity of the engine system to such changes, and a sophisticated numerical model, especially when 
used in conjunction with measured data, is a highly effective analytical tool. 

In the current application, a general-purpose CFD code named PHOENICS, developed by CHAM Inc., is 
used to model the temperatures, pressures, and velocities in the SSME’s High Pressure Fuel Turbopump (HPFTP) 
aft-platform seal cavity for a variety of boundary conditions and geometries. This cavity is located downstream of 
the fuel pump’s second turbine disk, between the disk and the aft platform seal (Figs. 1 to 4). It is an annular cavity 
where 1400°R combustion products and 150°R coolant hydrogen mix in a complex flow pattern and then are vented 
into the pump’s turbine exhaust. An understanding of the flow field in this cavity is critical since there are at least 
two known problems in the High Pressure Fuel Pump which may be linked to the environment in this region. 
Specifically, these problems are (1) cracking of the second stage turbine blade shanks, and (2) hot gas leakage into 
the stack behind the aft platform seal (Fig. 4). The first problem, blade cracking, can severely limit the time a pump 
can operate before it must be rebuilt. The second problem, that of hot gas leakage, is potentially more severe since, 
in the extreme, it may cause the pump to shut down prematurely if the temperatures or pressures in the coolant liner 
behind the aft-platform seal exceed certain redlines. 

Accordingly, the primary purpose of the present analysis is to investigate the two problem areas mentioned 
above. In doing so, the study addresses the following questions: 

1) How severe is the temperature gradient in the region where the turbine blades are cracking? 

2) What would be the temperature of any fluid which leaked from the cavity into the coolant liner? 

The analysis addresses these questions, not only for the pump operating under normal conditions, but also 
for a range of off-design conditions since even a slight departure from the norm might have a radical effect on the 
flow pattern and temperatures in the aft-platform seal cavity. As such, the broad objective of this study is to develop 
a model flexible enough that it can examine the effect that boundary parameters such as clearances, pressures, and 
flowrates have on the flow pattern and temperatures in the cavity. Such a model must be general enough that it can 
support future analytical and experimental investigations of the HPFTP aft-platform seal cavity. 
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Figure 2. SSME High Pressure Fuel Turbopump. 


Figure 3. HPFTP turbine flow paths. 



Figure 4. Aft-platform seal cavity. 






II. PROBLEM DESCRIPTION 


The region being modeled is the aft-platform seal cavity downstream of the second-stage turbine in the 
Space Shuttle’s HPFTP. A close-up of the aft-platform seal cavity is provided in Figure 4. General views of the 
j shuttle engine, the fuel pump, and the fuel pump turbine section are given in Figures 1 to 3. The dashed lines in the 
} close-up view, Figure 4, represent the limits of the problem as specified in the model. The fluid properties and either 
I the pressures or the flowrates at the boundaries must be input into the program. Unfortunately, the only available 
' measurements of these parameters (i.e., temperature and pressure) are far removed from the inlets and exits of the 
aft-platform seal cavity. As such, the boundary conditions chosen as inputs rely heavily on an existing one- 
dimensional analysis of the HPFTP and must be used with caution [1]. 

i Inspection of Figure 2, the HPFTP, will show that the aft-platform seal cavity is an axisymmetric annular 

cavity defined by stationary walls on one side and a rotating disk on the other. Flow enters the cavity through two 
inlets, one at the inner radius of the cavity and one near the outer radius of the disk. The flow leaves the region via 
i the gap between the outer radius of the aft-platform seal and the blade lips. At high rpm (up to 37,000) the flow is a 
i turbulent mixture of hydrogen and water at temperatures ranging from approximately 140°R to possibly as high as 
the turbine exhaust at 1700°R. Flowrates are on the order of 1 lbm/sec and pressures are in the range of 4000 psi. 

I 

The inlets and exits of the aft-platform seal cavity are described qualitatively below. The specific numbers 
used in this study, e.g., flowrates, pressures, etc., and the assumptions used in defining these numbers, can be 
found in the section on numerical model set-up. 

| 

J A. Inlets 

| 1 . Coolant Inlets 

I At the inner radius of the cavity, approximately 0.3 lbm/sec of liquid hydrogen flows into the aft-platform 

| seal cavity through a labyrinth seal. The source of this hydrogen is the coolant circuit which is fed by the discharge 
1 of the HPFTP (Fig. 3). In the two-dimensional model, this flowrate is calculated implicitly based on the pressure 
| drop through the labyrinth seal. In the three-dimensional model, in the interest of computational economy, the 
1 coolant flowrate through the labyrinth seal is not calculated internally, but is simply set to the value predicted by the 
two-dimensional model operating with the same average clearances and flowrate through the blade shanks. 

! 


2. Hot Gas Inlet at the Blade Shanks 

One wall of the aft-platform seal cavity is formed by the rotating disk upon which are mounted the second 
! stage turbine blades. At the periphery of this disk, a mixture of coolant hydrogen and combustion products enters the 
| cavity through the gap between the shank of one turbine blade and the next (Fig. 5). Since there are 58 blades in the 
j second stage disk, there are, accordingly, 58 holes available for this hot gas mixture to flow through into the aft 
platform seal cavity from the high pressure side of the turbine disk. The flow pattern of the fluid entering through 
these holes is complex since the shanks of the blades are curved and the disk itself is rotating at up to 37,000 rpm. 

In modeling this inlet, the 58 separate streams entering through the disk have been “smeared” in the circum- 
ferential direction into a single, continuous axisymmetric source. The flowrate and fluid properties at this inlet are 
prescribed based on predicted values, and the angular velocity of the fluid entering the cavity through these passages 
is assumed to have the same angular velocity as the disk. 
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HPFTP second stage turbine disk with blade: 
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| B. Exits 

j 

1 . Exit Gap Between the Outer Diameter of the Aft-Platform Seal and the Blades 

In this study, the single most important parameter which affects the flow pattern in the aft-platform seal 
cavity is the gap between the outer diameter of the aft platform seal and the lip of the turbine blades. This gap is very 
Ismail, on the order of one hundredth of an inch, and it supports a high pressure drop of over 500 psi between the 
j aft-platform seal cavity and the turbine exhaust. Any slight variation in this gap clearance will have a strong effect 
! on the total flow and overall flow pattern in the cavity. In general, the actual flow exiting through this gap at a given 
location will respond to changes in the overall turbine discharge pressure, the circumferential variation in turbine 
discharge pressure, and any changes in the width of the gap. The latter could be due to a number of different causes, 
including: sideloads, dynamics, machining tolerances, eccentricity, or thermal expansion. 

In the model, the exit pressure outside the gap is fixed at the best estimate for the turbine discharge pressure. 
The pressure drop across this exit is then related to the flowrate based on a loss coefficient times the local dynamic 
head. 

* 

i 

2. Secondary Exit Hole 

l 

In one of the test runs discussed in this report, the aft-platform seal is modeled with a second exit in order to 
simulate a postulated leak. The leak was assumed to be around the bolts which secure the aft-platform seal to the 
lift-off seal stack (Fig. 4). The hole size, loss coefficient, and exit pressure of this second exit were chosen such that 
the resulting calculated flowrate would be approximately 0.2 lbm/sec. The leak rate of 0.2 lbm/sec was chosen 
because it is the maximum flowrate which could be leaking past the bolts. This last conclusion is based on experi- 
mental measurements of the pressure drops in the coolant liner cavity which is downstream of the postulated leak . 


| III. NUMERICAL MODEL SET-UP 

I 

I 

CHAM Inc.’s general purpose computational fluid dynamics code, PHOENICS [2], has been employed for 
j all the numerical studies described herein. To use PHOENICS, special purpose “satellite” and “ground station” 
sub-programs must be formulated whereby the built-in features can either be turned on or off or modified, as neces- 
sary. One set of the sub-programs adapted specifically for the HPFTP aft-platform seal cavity three-dimensional 
studies is listed, in full, in Appendix A. Full listings of the other adapted sub-programs used in this study are given 
j in a separate CHAM report [3]. All of these sets of sub-programs are extensively annotated (via built-in 
| “COMMENT” statements) so as to make then self-explanatory when read in conjunction with the PHOENICS 
! User’s Manual [4]. Consequently, no detailed line-by-line description is given here; however, the most relevant 
features are described below. 

The two-dimensional calculations described herein have been performed by using the two-dimensional y/z, 
polar coordinate option of the code. Figure 6 shows the selected two-dimensional grid distribution. There are 1 120 
control cells, with 40 and 28 cells in the radial (IY) and (IZ) directions, respectively. Due to the (initially) assumed 
cyclic symmetry of the problem, only one control cell is required in the circumferential (IX) direction. However, to 
enable correct account to be taken of the wall shear stresses acting on the fluid entering between the blade shanks, 
the circumferential extent of the calculation domain is taken to be equal to the space between 2 consecutive blades 
(i.e., an angle of 1/58 X 2 it deg, where 58 = total number of blades). 

In the three-dimensional calculation, the full three-dimensional x/y/z coordinate capabilities of PHOENICS 
were employed. The identical y/z grid distribution of the 2-dimensional calculations was retained with, in addition, 
8 cells in the circumferential (IX) direction, such that a total of 8 X 40 X 28 = 8960 control cells is used. 
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As depicted in Figure 3, the “cold” liquid hydrogen coolant enters axially, through the labyrinth seal, at the 
inner radius of the cavity. This “cold” hydrogen then joins with the mixture of “hot” hydrogen and water that flows 
into the cavity from between the blade shanks located at the outer radius of the rotating disk. The combined streams 
of fluid then exit beneath the blade lips, as also shown in Figure 3. 


A. Assumptions/Model Details 

The major assumptions and salient features of the physical models and the boundary conditions employed 
are described below. 

1 ) All boundary surfaces (both stationary and rotating) have been assumed to be adiabatic . 

2) The hydrogen and water mixture are treated as a single homogeneous fluid with mixture properties 
(density and laminar viscosity) and temperature deduced from the calculated mixture enthalpy and specified 
hydrogen and water property curve fit data as described in Appendix B . 

3) The turbulence effects are presented by way of the two-equation (k-€) model of turbulence. In this model, 
two parameters, viz: the turbulence kinetic energy, k, and its dissipation rate, e, are computed from differential 
transport equations. Thus, it has the capability of representing both the local and history effects. The effective 
viscosity is expressed as: 


Feff = + p k 2 /e 
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|vhere (x^ is the laminar viscosity, is an empirical constant and p is the local mixture density. In addition, four 
ijther empirical constants are assigned the values as recommended in original publications [4], 


4) All boundary surfaces of irregular shape are accommodated in the present calculations by use of “cell 
porosities.” In this approach, each control cell is characterized by a set of fractions, in the range from 0 to 1 . These 
fractions determine the proportion of the cell volume which is available for flow from the cell to its neighbor in a 
*iven direction. This practice is much more rigorous and accurate than the practice of using rectangular steps. 

5) The wall shear stress is calculated by using the conventional wall functions which are based on the 
assumption of the logarithmic law of the wall . For partially blocked control cells, the wall stress is calculated for the 
projected surfaces parallel to the velocity components. 

It should be noted that the PHOENICS (1981 version) built-in process for determining wall shear stress 
is restricted to a finite number of special regions, to be set via the satellite subroutine. For the complex aft-platform 
seal geometry, many such special regions would be necessary, in excess of the built-in maximum, and a special 
PHOENICS user subroutine program was written for the current problem to overcome this restriction. This user 
isub-program (GWALL) performs the identical job as the built-in PHOENICS “WALL” subroutine but is used via 
the PHOENICS ground station. A listing of GWALL is included in Appendix A. 

! 

1 6) In PHOENICS, an iterative finite-difference solution procedure is employed to solve the governing 

differential equations together with the above mentioned relations. The method is based on a fully implicit, con- 
| servative formulation. As a result there is no restriction on the selection of the grid and the magnitude of the time 
steps. 

The variables calculated and/or solved for (and printed) in the seal cavity flow calculation include the 

following: 

1 a. The fluid velocities in the 3 coordinate directions 

b. The mixture enthalpy and deduced temperature 
| c. The (mass) concentration of water vapor 

d. The turbulent kinetic energy and its dissipation rate 

e. The static and total pressures 

f. The mixture density and separate densities of both the hydrogen and water 

g. The effective viscosity. 


7) Boundary conditions are: 


a. Prescribed mass flowrate, velocities, enthalpy, mixture ratio, and turbulence parameters at all inlets 
j except for the two-dimensional solutions, in which case the flowrate through the labyrinth seal was computed based 
on a prescribed inlet pressure. 


b. Prescribed exit pressure at all outlets, with the pressure drop related to the flowrate based on a speci- 
fied loss coefficient times the dynamic head. 
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c. The incoming fluid enclosed between the blade shanks is assumed to rotate at the same speed as the 


adjacent disk surface. 


8) The (phase change) freezing of the water is not accounted for; any water at temperatures below freezing is i 

given the properties (density, etc.) of liquid water at freezing. ! 

9) The effects of viscous heating have been ignored. J 

IV. TWO-DIMENSIONAL TEST RUNS ^ 


Three different two-dimensional test cases were run. The first of these was considered to be the basecase 
using the best estimate of the average conditions for the pump operating at the full power level (FPL). A second test 
run was made with a reduced amount of coolant entering through the labyrinth seal in order to determine the sensi- 
tivity of the solution to the ratio of hot gas flowing in at the blades relative to the hydrogen entering at the labyrinth 
seal. Finally, a third two-dimensional test run was made in order to see what effect a postulated leak through the 
stack bolts would have on the calculated cavity temperatures and flows. These three test runs and results are 
described in more detail below. 


A. Two-Dimensional Test Runs: Boundary Conditions 

1 . Basecase 2-D 

The basecase two-dimensional run uses boundary conditions and operating clearances taken from a one- 
dimensional flow analysis provided by Lockheed, Inc. [1]. These boundary conditions are tabulated in Table 1 . It 
should be noted that for this particular run, the boundary condition specified at the labyrinth seal is that of a pre- 
scribed pressure boundary from which the flowrate is then deduced based on the following relationship [5]: 


MASSFLOW = FC * AREA * SQRT( ( RH0*P0 ( 1- ( PN/PO ) **2 ) ) / 

(NUMBER OF TEETH + ALOG(PO/PN) ) ) 

WHERE PO = UPSTREAM PRESSURE; PN = DOWNSTREAM PRESSURE; FC = FLOW COEFF. 


(Note that for the basecase test run, the above equation when coupled with the PHOENICS two-dimensional model 
predicts a slightly lower flowrate through the labyrinth seal (0.26 lbm/sec versus 0.36 lbm/sec) as compared to the 
Lockheed one-dimensional model predictions.) 


2. Reduced Coolant (Labyrinth) Flow 

In the second run, the basecase two-dimensional model was modified by reducing the clearance at the outer 
diameter of the aft-platform seal while leaving all the other boundary conditions, including the hot gas flowrate, the 
same. When the gap size is reduced, the pressure in the cavity goes up and the coolant through the labyrinth seal 
decreases. The purpose here was to determine the effect that a reduction in coolant flow would have on the tempera- 
ture field in the cavity. 
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3. Leak Through the Stack Bolts 


The final two-dimensional run of the current study simulated a 0.2 Ibm/sec leak through the stack bolts. The 
boundary conditions for this run were the same as the basecase but with a “hole” at the location shown in Figure 4. 
The loss coefficient at this hole and the hole size were chosen such that they dictated a leak rate of approximately 0.2 
lbm/sec. 


TABLE 1. TWO-DIMENSIONAL BOUNDARY CONDITIONS 

Reduced 


Variable 

Basecase 

Coolant 

Leak 

Rotational speed of the disk (RPM) 

37,000 

37,000 

37,000 

Gap size at the labyrinth seal (in.) 

0.1069 

0.1069 

0.1069 

Total flow area (360°) between the blade shanks (in. 2 ) 

3.877 

3.877 

3.877 

Clearance between the aft-platform seal and blades (in.) 

0.0108 

0.0102 

0.0108 

Loss coefficient at the exit near the blade shanks 

1.5 

1.5 

1.5 

Enthalpy of the H 2 upstream of the labyrinth seal (Btu/lbm) 
(Resultant calculated temperature - degrees Rankine) 

278.3 

(145°R) 

278.3 

(145°R) 

278.3 

(145°R) 

Enthalpy of H 2 and H 2 0 entering through the blades (Btu/lbm) 
(Resultant calculated temperature - degrees Rankine) 

3558 

(1466°R) 

3558 

(1466°R) 

3558 

(1466°R) 

Density of the H 2 upstream of the labyrinth seal (lbm/ft 3 ) 

3.574 

3.574 

3.574 

Density of H 2 and H 2 0 entering through the blades (lbm/ft 3 ) 

0.931 

0.931 

0.931 

Mass flowrate of H 2 and H 2 0 entering past the blades (lbm/s) 

3.649 

3.649 

3.649 

Mass fraction of H 2 0 entering through the blades 

0.474 

0.474 

0.474 

Pressure at the turbine discharge (psi) 

3582 

3582 

3582 

Pressure at the labyrinth seal inlet (psi) 

4254 

4254 

4254 

Loss coefficient for the second (leak) exit 

1.5 

1.5 

1.5 

Total flow area at the second (leak) exit (in. 2 ) 

0 

0 

0.00019 
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B. Two-Dimensional Test Runs: Results and Observations | 

i 

1 . Basecase I 

According to the model, when the HPFTP is operating at full power, centrifugal force dominates the flow | 
field and pressure field in the aft-platform seal cavity. This is not surprising when one considers that at 37,000 rpm 
and a radius of 4.5 in., the hot gas exits the blade shanks with a centrifugal force equal to approximately 175,000 ■ 
g’s. Figures 7 and 8 show that there is virtually no penetration of the hot gas down into the aft-platform seal cavity 
and, as a result, the temperature in the cavity remains cold, at approximately 375°R (-85°F). The flow pattern in the \ 
main cavity consists of two large co-rotating vortices which maintain the cavity at a relatively uniform temperature. I 
The liquid hydrogen which enters through the labyrinth seal at the inner radius of the cavity flows radially outward | 
along the face of the disk and then abruptly merges with the hot fluid stream exiting from between the shanks. While I 
it appears from the drawing that the cold flow then recirculates, in fact all of the coolant which enters through the ! 
labyrinth seal must, by continuity, mix and then exit with the hot stream, resulting in the sharp temperature gradient | 
at the blade shanks especially evident in the close-up view provided in Figure 8. The actual local gradients that the 
blade shanks would see would be more severe than predicted here since, in the model, the hot gas flow is treated as 
an axisymmetric source which would tend to smooth out the temperature gradients at the trailing edge of the blade 
shanks. In reality, there are 58 blades shanks between which the hot gas flows into the cavity. The cold fluid which 1 
is slung off the disk, up behind the trailing edge of the blade shanks, will be sheltered from the hot flow entering 
from between the blade shanks. As a result, the local mixing of hot and cold fluid will be delayed, and the local { 
temperature gradient will be even more severe than shown here. 

i 

2. Reduced Coolant (Labyrinth) Flow I 

During the course of the study the question was raised as to what would happen if the proportion of coolant to 
hot gas flow were different that that predicted by the one-dimensional model used to define the boundary conditions 
[1]. In order to answer this question, the boundary conditions in the model are manipulated in a somewhat contrived ; 
manner in order to change the proportion of hot gas flow to coolant flow, viz: the coolant flowrate is reduced by 
slightly reducing the clearance between the aft-platform seal and the blade lips. The result is that a reduction of only 
six ten-thousandths of an inch (6 percent) of this clearance reduces the coolant flow by over half. This change, 
however, has little effect on the flow field in the cavity. As with the basecase, the flow in the cavity remains 
dominated by the centrifugal force. As shown in Figures 9 and 10, the temperature in the cavity has risen by only 
approximately 150 deg, up to 525°R, which is a moderate increase when compared with the hot gas inlet tempera- 
ture of 1466°R. The conclusion is that the flow field and temperature field of the aft-platform seal cavity is relatively 
insensitive to the amount of coolant entering at the labyrinth seal relative to the amount of hot gas mixture entering 
through the blade shanks. However, the pressure and coolant flowrate are extremely sensitive to the exit clearance at 
the outer diameter of the aft-platform seal for a fixed hot gas inlet flow. 


3. Leak Through the Stack Bolts 

For the third two-dimensional test run, a “hole” is simulated underneath the bolts which secure the lift-off 
stack. The rationale behind such a study is that a flow leaking past these bolts into the coolant liner might be one 
explanation for the erratic temperatures and pressures sometimes recorded in the coolant liner. The exit area and loss 
coefficient at this hole are adjusted so that the calculated leakage rate is 0.2 lbm/sec. The flowrate of 0.2 Ibm/sec 
comes from the best estimate of the upper limit of what the leak rate could be, based on the known temperature and 
pressure measurements in the liner [6]. 

Figures 1 1 and 12 show that a leak of 0.2 lbm/sec through the stack bolts does not dramatically change the 
flow field or temperature field as compared with the no-leak, baseline case. The temperature of the main cavity and 
the fluid leaking out past the bolts remains relatively unchanged at around 375°R (-85°F). 
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Figure 7. Two-dimensional basecase results 
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Figure 8. Two-dimensional basecase results, expanded view. 
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Figure 1 1 . Two-dimensional 0.2 lbm/s leak flow results 
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Figure 12. Two-dimensional 0.2 lbm/s leak flow, expanded view 









C. Convergence Characteristics and Computer Time 

Numerical solutions of flows involving rotating boundaries are notoriously slow to converge for a variety of 
reasons (not to be discussed here) and so it was deemed essential that careful checks be made to ensure that the 
PHOENICS solutions being obtained were meaningful. Thus, before the two-dimensional production runs 
described above were fully completed, a series of test calculations were performed to ensure that the solutions were 
converged to an acceptable degree. To this end, various runs were made for the basecase setup with different initial 
guess/starting solutions that were quite extensive. The results of these investigations are presented and discussed in 
Appendix C. As shown in the latter, the PHOENICS solutions are clearly converging to an identical solution in each 
case, as should (and must) be expected. 

As depicted in Appendix C, the two-dimensional basecase was run for a total of 500 sweeps, at which time 
all calculated monitor flow variables had settled to an acceptable degree (Figs. C-l to C-8). All the other 2- 
dimensional runs reported here were restarted from this basecase solution (i.e., the initial fields for the starting of 
the iterative calculation procedure were taken to be the basecase solution, rather than some simple initial guess) and 
then run on until, again, the solution monitor values were suitably settled. This usually required another 150 to 200 
sweeps, at most. Computer times for these restart runs were approximately 35 CPU minutes on CHAM’s Perkin 
Elmer 3251 mini-computer. 

All the 3-dimensional calculations described in the next section were also restarted from the 2-dimensional 
basecase solution which was symmetrically duplicated in the circumferential IX-direction. Again, converged 
solutions then took approximately 150 to 200 more sweeps and required approximately 5 CPU hours of computer 
time on the Perkin-Elmer 3251 machine. 


V. THREE-DIMENSIONAL TEST RUNS 


The disadvantage of the preceding axisymmetric analysis is that, by definition, it does not include the three- 
dimensional effects either known or suspected to exist in the pump. One of the most important of these asymmetries 
left unaccounted for by the two-dimensional analysis is the circumferential variation in pressure which has been 
measured downstream of the exit of the fuel turbine. This exit pressure serves as one of the boundary pressures 
which regulates the flow in the aft-platform seal cavity. In addition to this known pressure variation, there may be 
variations in clearances or other parameters which could radically alter the flow pattern in the cavity. As such, a 
three-dimensional model is an essential tool for a proper study of this cavity. As a starting point, three different 
three-dimensional cases were run and are presented here. The first is the basecase which uses the same set of 
flowrates, fluid properties, and clearances as used in the two-dimensional basecase. The only difference between 
the two is that the three-dimensional basecase also includes a prescribed asymmetrical turbine exit pressure based on 
pressure measurements taken during a full scale test of the shuttle engine. The second three-dimensional case was 
set-up to simulate a 0.003 in. shift in the rotor position with a corresponding change in the clearance at the labyrinth 
seal and at the exit gap between the aft-platform seal and the blade lip. This shift is relative to the average labyrinth 
seal clearance of 0.003 in. and the average exit gap of 0.0108 in. The last three-dimensional run presented here 
simulates a relatively large eccentricity of the aft-platform seal alone, such that the exit clearance is skewed to one 
side by 0.0081 in., which is 75 percent of its average clearance. 

A. Three-Dimensional Test Runs: Boundary Conditions 

1 . Basecase (Geometrically Axisymmetric with Asymmetric Exit Pressures) 

As its boundary conditions, the basecase three-dimensional run uses the same operating clearances, 
flowrates, pressures, mixture ratios, and enthalpies, etc., as used by the two-dimensional basecase analysis. The 
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only exception is that the exit pressure of the turbine is no longer uniform but varies circumferentially based on data 
taken during Rocketdyne’s engine test 902-279 [7], These boundary conditions are the best estimate of the operating 
conditions in the fuel pump at full power (109 percent). The specific numbers used for this run, and for the subse- 
quent three-dimensional runs, are listed in Table 2. 


2. Eccentric Rotor (Rotor Shift of 0.003 in.) 


The Eccentric Rotor (0.003 in.) case was set up to simulate the effect that a rotor shift of 0.003 in. would 
have on the flow field in the cavity. The shift of 0.003 in. was chosen because it is an upper limit on the distance the 
rotor can shift before the shaft starts rubbing against the labyrinth seal. Such a rotor shift in a given direction would 
open up the exit clearance between the aft-platform seal and the blade lips, while at the same time it would close 
down the clearance at the labyrinth seal. This effect is simulated in the model by, on the one hand, directly adjusting 
the clearances at the outer diameter of the aft-platform seal and, on the other, by adjusting the flow rate at the 
labyrinth seal. All the other inputs remain the same as for the three-dimensional basecase. 

TABLE 2. THREE-DIMENSIONAL BOUNDARY CONDITIONS 


Variable 


Rotational speed of the disk (RPM) 

Flowrate at the labyrinth seal (lbm/sec) 1:00 

2:30 

4:00 

5:30 

7:00 

8:30 

10:00 

11:30 

Total Mass Flowrate 

Total How area (360°) between the blade shanks (in. 2 ) 

Clearance between the aft-platform seal and blades (in.) 1:00 

2:30 

4:00 

5:30 

7:00 

8:30 

10:00 

11:30 

Total Area 

Loss coefficient at the exit near the blade shanks 

Enthalpy of the H 2 entering at the labyrinth seal (Btu/lbm) 
(Resultant calculated temperature - degrees Rankine) 

Enthalpy of H 2 and H 2 0 entering through the blades (Btu/lbm) 
(Resultant calculated temperature - degrees Rankine) 

Density of the H 2 entering at the labyrinth seal (lbm/ft 3 ) 
Density of H 2 and H z O entering through the blades (lbm/ft 3 ) 
Mass flowrate of H 2 and H 2 0 entering past the blades (lbm/s) 
Mass fraction of H 2 0 entering through the blades 

Pressure at the turbine discharge (psi) 1 :00 

2:30 

4:00 

5:30 

7:00 

8:30 

10:00 

11:30 

Average Exit Pressure 


Basecase 

Rotor 

Eccentricity 
- 0.003 in. 

A ft -PI at form 
Eccentricity 
= 0.0081 in. 

37,000 

37,000 

37,000 

0.0323 

0.0095 

0.0323 

0.0323 

0.0323 

0.0323 

0.0323 

0.0551 

0.0323 

0.0323 

0.0646 

0.0323 

0.0323 

0.0551 

0.0323 

0.0323 

0.0323 

0.0323 

0.0323 

0.0095 

0.0323 

0.0323 

0.0000 

0.0323 

0.258 

0.258 

0.258 

3.877 

3.877 

, 3.877 

0.0108 

0.0129 

0.0165 

0.0108 

0.0108 

0.0108 

0.0108 

0.0087 

0.0051 

0.0108 

0.0078 

0.0027 

0.0108 

0.0087 

0.0051 

0.0108 

0.0108 

0.0108 

0.0108 

0.0129 

0.0165 

0.0108 

0.0138 

0.1800 

0.307 

0.307 

0.307 

1.5 

1.5 

1.5 

278.3 

278.3 

278.3 

(145°R) 

(145°R) 

(145°R) 

3380 

3380 

2280 

(1466°R) 

(1466°R) 

(1466°R) 

3.574 

3.574 

3.574 

0.931 

0.931 

0.931 

3.649 

3.649 

3.649 

0.474 

0.474 

0.474 

3451 

3451 

3451 

3541 

3541 

3541 

3697 

3697 

3697 

3622 

3622 

3622 

3606 

3606 

3606 

3592 

3592 

3592 

3476 

3476 

3476 

3481 

3481 

3481 

3558 

3558 

3558 
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3. Eccentric Aft-Platform Seal (Aft- Platform Seal Shift of 0.0081 in.) 


The third three-dimensional test run models the flow field for a highly eccentric (75 percent) aft-platform 
seal. In this run the clearance at the gap between the aft-platform seal and the blade lips is adjusted so that it models 
what the gap would be if the aft-platform seal had moved laterally 0.0081 in. in the 1 1 :30 direction. Note that the 
rotor itself has not moved but is still concentric with the labyrinth seal so that the gap between the labyrinth seal and 
the rotor axle remains at a uniform 0.003 in. (In general, the clocking positions used in this report correspond to the 
convention adopted by Rocketdyne in Reference 7, however, in this particular test run the decision to move the 
aft-platform seal in the 1 1 :30 direction is arbitrary, and is based on convenience rather than any physical justifica- 
tion.) The choice of the magnitude of the eccentricity is also somewhat arbitrary but the reasoning behind the shift of 

0.0081 in. was the desire to choose a large aft-platform eccentricity in order to observe extreme effects. An aft- 
platform shift of 0.0081 in. is 75 percent of the total aft-platform seal clearance. 


B. Three-Dimensional Test Runs: Results and Observations 


1. Basecase 

A comparison of the three-dimensional basecase results (Figs. 13 to 21) with the two-dimensional basecase 
results (Figs. 7 and 8) shows that the addition of an asymmetric pressure distribution at the exit of the turbine has had 
little effect on the flow pattern in the aft-platform seal cavity. While some evidence of the influence of the external 
pressure distribution can be seen at the outer diameter of the disk near the blade shanks (e.g. , Fig. 1 6), this effect is 
small; toward the center of the cavity the results are nearly identical to the two-dimensional solution. At the 
flowrates and small clearances of the aft-platform seal cavity running at full power, a circumferential pressure 
difference of 220 psi as modeled here represents only a fraction of the over 600 psi pressure drop between the 
aft-platform seal cavity and the turbine exhaust. As a result, the 220 psi circumferential variation on the outside of 
the cavity has little effect on the flow pattern inside. In addition, even with the circumferential variation in turbine 
exhaust pressure, the centrifugal force in the aft-platform seal cavity still dominates the flow such that the influence 
that is felt due to the pressure variation is confined to the periphery of the cavity. For an example of this effect, 
examine the lines of constant temperature given in the close-up view in Figure 17. 


2. Eccentric Rotor (Rotor Shift of 0.003 in.) 

Perhaps the most notable feature of the aft-platform seal cavity flow field (Figs. 22 to 30) with a 0.003 in. 
eccentric rotor is the small change as compared to the three-dimensional basecase with its centered rotor. Even with 
an eccentric rotor, the temperatures in the cavity have risen just 75°R, indicating only a slight increase in the heat 
transferred down into the cavity. Again, the only significant effect is felt at the outer diameter of the turbine disk 
where, at the 5:30 clock position, the hot gas actually flows down into the cavity causing a local hot spot. This hot 
spot will be felt by the blade shanks once per revolution, with a corresponding cooling in between. In general, 
therefore, a rotor shift of 0.003 in. results in a slight warming of the average cavity temperature, and a cyclical 
variation of temperature at the outer diameter of the disk of approximately 600°R. 


3. Eccentric Aft-Platform Seal (Aft-Platform Seal Shift of 0.0081 in.) 

Of the six different two-dimensional and three-dimensional test runs investigated during this study, the most 
dramatic results (Figs. 31 to 39) come from running the model with an aft-platform seal that has been shifted to one 
side by 3/4 of the exit clearance (i.e., by 0.0081 in.). With the exit gap substantially closed down on one side, the 
hot gas which would normally exit through that gap must, instead, exit at a different location. The centrifugal force 
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Figure 18. Three-dimensional basecase results: temperature (end view) 
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Figure 27. Three-dimensional eccentric (0.003 in.) rotor: temperature (end view) 
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Figure 28. Three-dimensional eccentric (0.003 in.) rotor: Mass concentration 
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Figure 32. Three-dimensional eccentric (0.0081 in.) aft-platform seal: vectors (close-up). 
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Figure 33. Three-dimensional eccentric (0.0081 in.) aft-platform seal: vectors (end view). 
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Figure 34. Three-dimensional eccentric (0.0081 in.) aft-platform seal: temperature 
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Figure 35. Three-dimensional eccentric (0.0081 in.) aft-platform seal: temperature (close-up). 



Figure 36. Three-dimensional eccentric (0.0081 in.) aft-platform seal: temperature (end view) 
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Figure 38. Three-dimensional eccentric (0.0081 in.) aft-platform seal: static pressure 






























on this hot gas still works to confine it to the outer radius of the cavity. However, the pressure differences in the 
radial direction are, in this case, becoming large enough that they are forcing more and more hot gas down into the 
I cavity. This is clearly evident in the velocity diagrams (Figs. 31 and 32) where, at the 5:30 clock position, there is a 
strong inward flow of hot gas down into the cavity. The temperature profiles also indicate a dramatic increase in hot 
gas in the cavity. The temperatures at the center of the cavity are now up to 675°R (2 1 5°) which is 300° warmer than 
for the basecase. As with the other three-dimensional runs, the most pronounced effects can still be seen at the outer 
j radius of the cavity. Here at the outer radius of the disk near the blade shanks, there is a large circumferential 
j variation in both temperature and pressure with the temperature cycling 600°R and the pressure varying by 20 psi. 

One further observation on the results from this test run has to do with the static pressure. The observation is 
that for a 0.0081 in. eccentric aft-platform seal, the pressure in the cavity has gone up by 80 psi relative to the 
! three-dimensional basecase. The implications of this pressure rise are not simple to determine. The difficulty lies in 
the fact that such a pressure in the cavity would reduce the flow rate through the labyrinth seal, which, for the 
three-dimensional model is (numerically) fixed based on the flowrates calculated earlier in the two-dimensional 
basecase. Since the total exit areas and average turbine discharge pressure for all the three-dimensional runs are the 
same as for the two-dimensional basecase, this is a reasonable assumption. In this final three-dimensional case, 
however, the assumption leads to a contradiction. For the eccentric aft-platform seal case, using the flowrates from 
the two-dimensional basecase, the pressure at the exit of the labyrinth seal is calculated as being higher than the 
pressure at the labyrinth inlet. In other words, if in the three-dimensional model this boundary had been specified as 
l a fixed pressure instead of a fixed flowrate, the model would have predicted reverse flow through the labyrinth seal. 

1 That there could actually be reverse flow through the labyrinth seal is considered extremely unlikely. What would 
more likely occur is that an eccentricity in the aft-platform seal would raise the pressure in the aft-platform seal 
cavity, reducing both the hot gas flow past the blade shanks and the labyrinth seal flow. 

VI. SUMMARY OF THE CURRENT TEST RUN RESULTS AND OBSERVATIONS 


, The axisy mmetric computer model of the aft-platform seal cavity indicates that at 37 ,000 rpm the flow in the 

I aft-platform seal cavity is dominated by the centrifugal force caused by the rotating turbine disk. The disk drives a 
recirculating flow in the central region of the cavity, creating a core of nearly uniform temperature. In general, the 
1 temperature field throughout the aft-platform seal cavity is dictated primarily by convection (as opposed to conduc- 
! tion) as indicated by the fact that little heat from the hot gas at the periphery of the cavity is conducted down into the 
isothermal core. As a result, the core stays relatively cold, even when the coolant flowrate is reduced by over 50 
percent. 

: The most severe temperature gradient in the aft-platform seal cavity occurs at the outer diameter of the 

turbine disk, near the blade shanks. At this location, the hot gas entering from between the blade shanks mixes with 
j the coolant flow that is being slung off the face of the disk. The temperature difference between the two streams is 
over 1000°R. 

I The three-dimensional computer model of the aft-platform seal cavity shows that, for normal clearances and 

operating conditions, the flow field in the cavity is relatively insensitive to the circumferential pressure variation 
known to exist in the turbine discharge. The flow field is shown to be sensitive, however, to eccentricities of the exit 
^ gap between the aft-platform seal and the blade shanks. But in both cases, it is the centrifugal force which still 
dominates the flow pattern, such that any perturbation of the flow field or temperature field which results from either 
pressure changes or geometrical changes are, for the most part, confined by centrifugal force to the outer diameter of 
the cavity. 


49 



In addition to the above, the study also reveals that, for fixed flow through the blade shanks, the labyrinth 
flowrate is extremely sensitive to the exit area at the outer diameter of the aft-platform seal. While this result is 
somewhat misleading, since it is based on the unrealistic boundary condition of a fixed flowrate through the blade 
shanks, it nevertheless merits further consideration especially with regard to transient phenomena. Finally, as a 
related observation, the flowrate through the labyrinth seal is also sensitive to the eccentricity of the aft-platform 
seal clearance, even for a constant exit area. This sensitivity is something which has yet to be included in the current 
one-dimensional models of the flow through the pump’s turbine section. 


VII. CONCLUSIONS 


The results of the study summarized above provide the following insight into the specific problems which 
initiated the study, i.e., (1) the cracking of the HPFTP second stage blades, and (2) the suspected hot gas leakage 
into the coolant cavity behind the aft-platform seal bolts. 

As far as the blade cracking is concerned, the model has shown that the second stage blade shanks are 
subjected to varying degrees of thermal stress, both steady state and once per revolution. The severity of this 
gradient has been shown to be sensitive to asymmetries in the external pressure and to variations in the geometrical 
clearances. At the time of this writing, however, it is believed that the primary cause of cracking is not due to 
thermal effects but is the result of a very high mean mechanical stress coupled with the moderate thermal stress. The 
proposed solution to alleviate the cracking is to recontour the shank in the high stress area, to shot-peen the surface 
to reduce the surface mean operating stress, and to coat the shanks to reduce the thermal stress [8]. 

As for the variations in coolant liner pressure and temperature thought to be indicative of a leak into the 
coolant liner, they remain an enigma. In order to gain a clearer understanding of this problem, the fluid temperatures 
calculated in the current study will be used as an input to the thermal stress analysis of the hardware. Prior to this 
study it was believed that the temperatures in the cavity were on the order of 900°R hotter than predicted here [8], 
With a better estimate of the fluid temperature, the thermal stress analysis will be better able to predict the deforma- 
tion of the aft-platform seal and the other components neighboring the aft-platform seal cavity. This will, in turn, 
generate improved estimates of the clearances and flowrates in the region. 

The new flowrates estimated from the above will be fed back into the PHOENICS model for an improved 
analysis of the flow and temperature field in the cavity. Other changes which could be incorporated into the model 
would be to include the effect of heat transfer into the cavity and the viscous heating of the fluid itself, both of which 
will result in increases in the cavity temperature. 

In addition to further analytical studies and improvements, there are plans to build a fuel pump that has 
pressure and temperature measurements built into the aft-platform seal, the labyrinth seal, the lift-off seal stack, and 
the coolant liner [8]. The test data from this instrumented pump, in conjunction with the computer model predic- 
tions should greatly increase the level of understanding of the operating environment of the high pressure fuel pump 
aft-platform seal cavity. 
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APPENDIX A: PHOENICS COMPUTER CODE 
SATELLITE AND GROUND ADAPTATIONS 
FOR THE SPACE SHUTTLE MAIN ENGINE HPFTP 
AFT-PLATFORM SEAL CAVITY 3-D MODEL 
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CALL CR( PC , O . O , 1 . NX , 6, 6. 1, 2.NX.NY.NZ) 
CALL CR(PN.O.O. 1 .NX. 6, 6. 1. 2.NX.NY.NZ) 
CALL CR(PE.O.O. 1.NX, 6. 6. 1. 2 , NX , NY , NZ ) 
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960 CALL CR( PH , RAT , 1 , NX ,32,32, 20 . 28 , NX , NY , NZ ) 

96 1 C 

962 C ROW 33 

963 CALL CR( PC ,0.0, 1 , NX ,33,33, 1 , 16 . NX , NY , NZ ) 
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1002 CALL GET ( JH 1 , GH 1 . NY , NX ) 
CALL GET< JT 1 , GT 1 .NY , NX ) 
CALL GET(JCI.GCI.NY.NX) 
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APPENDIX B: PROPERTY CURVE FITS 


The individual enthalpy curves for water and hydrogen have been combined in order to calculate a mixture 
enthalpy, Enthalpy mix , defined as: 

Enthalpy mjx (T) = (Mass Ratio H 2 0) * Enthalpy Water (T) + 

(1-Mass Ratio H 2 0) * Enthalpy Hydrogen (T) 

I 

| This combined property curve is needed to be able to calculate the temperature of any given mixture of water 

and hydrogen in the aft-platform seal cavity, based on the mixture ratio and enthalpy calculated by the model. From 
the temperature are then calculated other fluid properties, such as density and viscosity. The curve fits used to 
compute these properties are depicted in Figures B-l to B-6. 


i 
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ENTHALPY (BTU/LBM) 


ENTHALPY OF WATER 1 


CURVE FIT I H (Btu/lbm) « -A2A. 5938 + .82AIAT + I . 3067x1 O'^T 2 

(A92<T<975R) 


CURVE FIT II H = 2289.552 - A. 577089T + 2 . 8l 52A9x I 0* 3 T 2 

(975R<T<ll8A.6R) 


CURVE FIT III H « -7363.69 + 6.913T 

(1 18A. 6R<J< 1 223- 3R) 


CURVE FIT IV H = 599-5881 - I.27177T + I . 369267x 1 0 _3 T 2 

(1223.3R£T<1281.AR) 


CURVE FIT V H = -307-5AA9 + I.I90721T 

(I28I.AR<T<I AOOR) 


STANDARD ERROR = A. 08 Btu/lbm 


These curves were fit to data taken from Thermodynamic Properties of 
S team , Joseph Keenan and Frederick Keyes, (New York: Wiley and Sons, 1 936) 
pp. 72-75. 



Figure B-l. 
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ENTHALPY (BTU/LBM) 


I 


i 


ENTHALPY OF HYDROGEN ** 


CURVE FIT I 

(170R<T<508R) H (Btu/lbm) = -5 . 92706x lO’^T 2 + 4.4 68995T 

-357.6903 


CURVE FIT II H = - 7 . 15694xlO'^T 2 + 3.557702T - 45.88906 

(508R<T<2000R) 


STANDARD ERROR = 4.39 Btu/lbm 


These curves are fit to data taken from the Hydrogen Technological 
Survey - Thermophysical Properties . Robert D. McCarty, (Washi ngton , D.C . : 
NASA Scientific and Technical Information Office, 1975) p. 472. 



Figure B-4. 
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VISCOSITY (LBM/FT-SEC x 10 5 ) 


VISCOSITY OF HYDROGEN 


Vise. ( Ibm/f t-sec xlO 5 ) = .4989 - 5.4575xlO~ 5 T 

+ 5.I824xIO'7 T 2 - 1.4948x10-'°t3 


STANDARD ERROR (xlo5) = 0.00047 Ibm/ft-sec 


This curve is fit to data taken from the Hydrogen Technological Surve 
Thermophys ical Properties , Robert D. McCarty, (Washington, D.C.: Scientific 
and Technical Information Office, NASA, 1975) p. 473. 


15 1 VISCOSITY OF 

HYDROGEN 

1.4 - 
1,3 - 
1.2 - 
1.1 - 









APPENDIX C: CONVERGENCE CHARACTERISTICS 


i 

I 


The insensitivity of the model to the initial values chosen for temperature and velocity is demonstrated by the 
solution sets given below. There were four different test cases run using identical boundary conditions but with the 
different guesses of velocity and temperature listed in the following table. Cases 1 to 3 were run to test the sensitivity 
of the solution to the initial temperature guess, and case 4 was run to check the sensitivity of the solution to the initial 
I choice of the velocity field. 

t 

' INITIAL FIELD VALUES 


j 

CASE 1 
CASE 2 
CASE 3 
CASE 4 


Temperature Theta Velocity 

Hot Guess Omega = 0.4 Disk Omega 

Best Guess Same as above 

Cold Guess Same as above 

Best Guess Omega = 0 

(Omega in radians/sec) 


After 500 sweeps, the values of velocity, temperature, and pressure, at a reference point in the middle of the 
cavity, have converged to the extent shown in Figures C-l to C-8. By 500 sweeps the constant pressure lines (Fig. 
C-2) for all four cases are very similar, as are the streamlines (Fig. C-8), and to a lesser extent the temperature 
profiles (Figs. C-l and C-l). Of the three, temperature is the slowest to recover from a poor initial guess. However, 
even with an initial temperature estimate 1000°R off the final values, the temperature at the monitoring point has 
converged to within 200° of the final value after 200 sweeps, and to within 75° of the final value after 400 sweeps. 
This is an acceptable convergence rate for our application, especially since the magnitude of the error is readily 
apparent from the slope of the temperature convergence curve (Fig. C-l). In the event that greater accuracy were 
required, the solution could be bracketed or else extended the necessary number of sweeps. 
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Figure C-7. Velocity field and streamlines, cases 1 to 4. 


132 


1 . REPORT NO. 

NASA TP-2685 


2. GOVERNMENT ACCESSION NO. 


3. RECIPIENT'S CATALOG NO. 


. TITLE ANO SUBTITLE 

Space Shuttle Main Engine High Pressure Fuel Pump 
Aft Platform Seal Cavity Flow Analysis 


5. REPORT DATE 

January 1987 


6. PERFORMING ORGANIZATION COO 


7. AUTHOR(S) 

S. A. Lowry and L. W. Keeton* 


8. PERFORMING ORGANIZATION REPOI 


9. PERFORMING ORGANIZATION NAME AND ADDRESS 

George C. Marshall Space Flight Center 
Marshall Space Flight Center, Alabama 35812 


10. WORK UNIT NO. 

M-544 


1 1. CONTRACT OR GRANT NO. 


12. SPONSORING AGENCY NAME AND ADDRESS 

National Aeronautics and Space Administration 
Washington, D.C. 20546 


13. TYPE OF REPORT & PERIOD COVE 

Technical Paper 


14. SPONSORING AGENCY CODE 


15. SUPPLEMENTARY NOTES 

Prepared by Structures and Propulsion Laboratory, Science and Engineering Directorate. 
*L. W. Keeton is employed by CHAM (NA) Inc., Huntsville, Alabama 


16. ABSTRACT 


A general purpose, three-dimensional computational fluid dynamics code named PHOENICS, developed 
by CHAM Inc. , is used to model the flow in the aft-platform seal cavity in the high pressure fuel pump of the space 
shuttle main engine. The model is used to predict the temperatures, velocities, and pressures in the cavity for six 
different sets of boundary conditions. The results are presented as input for further analysis of two known problems 
in the region, specifically: erratic pressures and temperatures in the adjacent coolant liner cavity and cracks in the 
blade shanks near the outer diameter of the aft-platform seal. 


17. KEY WORDS 

Space Shuttle Main Engine 
High Pressure Fuel Turbopump 
Aft Platform Seal Cavity 
Computational Fluid Dynamics 
Turbine Blades 


18. DISTRIBUTION STATEMENT 


Unclassified — Unlimited 



Subject Category 34 

19. SECURITY CL ASSIF. (of this report) 

20. SECURITY CLASS! F. (of this page) 

21. NO. OF PAGES 

22. PRICE 

Unclassified 

Unclassified 

137 

AO 7 


For gale by National Technical Information Service. Springfield, Virginia 221 

NASA-Langley, 19 











